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Scintillating crystal detectors offer potential advantages in low-energy low-background 



O ' experiments for particle physics and astrophysics. The GSO crystal is an interesting detec- 

O ' tor to explore for future neutrino physics experiments. The contributions to background 

>^ ' due to the various channels of intrinsic radio-isotopes from the ^^^Th and ^^^U series 

I , are identified and studied with time-correlation analysis and detailed fits to the spectral 

O , shape. Good agreement is achieved between measured and simulated spectra, indicating 

background suppression factors to the 10~^ — 10~^ level are possible. The procedures can 



•^ ■ be adopted for background understanding and suppression in other low-count-rate exper- 

^ ' iments where the dominant source of background is from internal radioactivity. Based 

on 1656 hours of data taking, limits on the double beta decay half-life for the various 
channels in ^^°Gd are derived. The limits for the neutrinoless and the Majoron modes are 
rj.oup(3 ^ 2.0(0.8) X 10^0 y and Tf^ > 3.0(1.9) x lO^^ y, respectively, at 68(90)% confi- 

2 2 

dence level. 
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1 Introduction 

Crystal scintillators have been widely used as radiation detectors as well as in medical and 
security imaging |l|]. Large detector systems, of scale several tens of tons in mass, have 
been constructed and made operational as electromagnetic calorimeters in high energy 
physics experiments 0. It is one of the few detector technologies (the other being loaded 
liquid scintillator) where a large range of nuclei can be turned into massive detector 
systems. 

The merits of scintillating crystal detectors in non-accelerator particle physics experi- 
ments have been recently discussed 0. Primarily, the high-Z composition of most crystals 
allows a compact design and provides large suppression of background due to ambient ra- 
dioactivity if a three dimensional fiducial volume definition can be realized. Experiment 
based on 100 kg of Nal(Tl) has been built for Dark Matter searches, producing some of 
the most sensitive results ||^. Another experiment toward 500 kg of CsI(Tl) for studying 
low-energy neutrino interactions using reactor neutrinos as probe is now under construc- 
tion 0. 

An interesting isotope for low-energy neutrino physics experiment is ^^°Gd, which has 
an isotopic abundance of 21.9% in natural gadolinium. It is a double beta decay can- 
didate [||. The lepton- number- violating neutrinoless double beta decay {0^(3(3) channel, 
which is sensitive to the Majorana neutrino mass terms, proceeds as: 

i60q^ ^ 160 j3y + 2 e" ( Q = 1.73 MeV ) . 

In addition, this isotope has recently been proposed as an attractive isotope for the detec- 
tion of low-energy solar neutrinos by providing a time-delay signature for background 
suppression and for tagging the flavor-specific z/cN-CC reactions [§]: 

Ue + ^'^"Gd ^ e- + 7's(prompt) + ^^^Tb* (Q = 169 keV) , 

followed by, in the 100 ns time-scale, 

leorj^* ^ leorj.^ ^ ^ ^^^ _ go ns . g^ = 63,7 ^eV) . 

The scintillating crystal GSO (Gd2Si05(Ce)), where gadolinium is 74.4% by weight, 
was first developed in the 1980's p. It is an attractive technology to realize an ac- 
tive detector based on gadolinium. Various studies have been performed at the proto- 
type [0 and detector |]1T| levels for calorimeter systems, as well as for double beta decay 



searches ||T2|, [13], Q and feasibility studies of solar neutrino experiments [T^ . 



The properties of GSO, together with a few commonly used crystal scintillators, is 
listed in Table |l|. It can be seen that besides neutrino physics, GSO may offer an attractive 



alternative to BGO in high energy physics and medical applications, having almost the 
same radiation length, but higher light yield, faster decay constant, and excellent radiation 
resistance. 

In this article, we report on the measurement and analysis of the major components 
of intrinsic radio-impurities in a commercially available GSO crystal^. These results are 
important to assess the feasibilities of GSO crystals in the low-event-rate experiments 
mentioned, and can serve as basis of future R&D efforts. In addition, the procedures 
can be adopted for background understanding and suppression in experiments dominated 
by internal radioactivity background. Lifetime limits on the various double beta decay 
channels for ^^°Gd are derived. 



2 Experimental Setup 

The schematic diagram of the experimental setup is depicted in Figure [l|a. The GSO 
crystal under study has a dimension of 3.8 x 3.8 x 18 cm^, weighing 1.744 kg and was 
wrapped with white Teflon tapes to improve the light collection. An outer layer of black 
vinyl tapes provided a light-proof layer. 

A 2-inch diameter photo-multiplier(PMT) with potassium-free glass was coupled to 
one of the end surfaces of the crystal. The whole crystal-PMT system was inserted into 
a well-shaped Nal(Tl) crystal scintillatoiQ with a mass of 14 kg, and viewed via a quartz 
window by a 5-inch photo-multiplier tube with potassium-free glass. The Nal(Tl) scin- 
tillator served as an active shielding device, vetoing background /?'s and 7's from the 
ambient environment, as well as the cosmic muons that penetrated the crystal detector. 
Surrounding the Nal-GSO-PMT active components was lead shieldings of 10 to 15 cm 
in thickness. The measurement was performed in a ground-level laboratory of a seven- 
storey building. The estimated overburden is about 12 meter- water-equivalent, sufficient 
to suppress the hadronic components of the cosmic-rays. A 40 x 50 x 1 cm^ scintillation 
plastic was placed on top of the detector outside the lead bricks for cosmic-ray muon ve- 
toes. As shown in Section |^, internal radioactivity becomes the dominant contribution to 
the background spectrum under this shielding configuration. Therefore, it is not necessary 
to make special additional efforts to suppress the neutron and radon induced background 
in the present work. 

The schematic diagram of the electronics and data acquisition systems is displayed in 
Figure |l|b. The CAMAC system was adopted. The energy and time for each event were 

•^Manufacturer: Hitachi Co. Ltd., Japan. 
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recorded. The energy measurement was made by using an analogue-to-digital converter^ 
which integrated the anode current of the PMT for a gate-time of 600 ns, while the 
timing measurement was done by a scaler^ counting on a 1 MHz clock. The scaler was 
read and then reset when the LAM (Look-At-Me) line was set active by the ADC. The 
timing accuracy was checked to be better than the time bin of 1 fis. Data were acquired 
and transferred through the CAMAC bus via a fast CAMAC Controller f\ to an interface 
processor running on the Linux operating system with a PC computer. The rate at a 
detector threshold of 50 keV after active vetoes and passive shieldings was 6.46 Hz. The 
dead-time of the system, which was mainly due to the conversion time of the CAMAC 
ADC (about 106 /is per conversion), is less than 0.1% 



3 Results and Analysis 

Events from the GSO crystal not in coincidence with the Nal(Tl) anti-Compton detector 
and the cosmic-ray veto panel were selected. The energy spectrum based on 1656 hours 
of data taking is shown in Figure |[ It can be seen that there are three distinct regions in 
the energy spectrum due to intrinsic radioactivity. The sharp peak at 0.394 MeV is due to 
the a-decay of the isotope ^^^Gd in the crystal. The nominal isotopic abundance of ^^^Gd 
is 0.2%, decaying by a-emission at 2.14 MeV and with a half life of 1.08 x 10^^ years [|16]. 
Events in the broad peak around 1 MeV are due to the a-decays of the radio-nuclides in 
the ^^^Th and ^ss-y ggj^jgg j^ the GSO crystal. These isotopes and their decay-daughters 
are primordial in origin and get into the crystal mostly through the contaminations in 
the raw materials. The high energy shoulder above 1.4 MeV and up to 2.0 MeV is due 
to the /3-emissions together with the subsequent prompt 7-rays from the radio-nuclides 
of these cascades. The bump at the 2.6 MeV region is due to the ^os^l 2.61 MeV 7-rays 
from ambient radioactivity. 

There is no notice- able peak at 1.46 MeV due to the 7-rays from the commonly found 
^°K, since a PMT with potassium-free glass was used. In the first trial measurement using 
another PMT with normal glass, a peak at this energy did show up prominently. 

The ^^^Gd a-peak at 394 keV served as a convenient reference standard to monitor 
and to correct for the gain instabilities of the measurement system in the course of data 
taking. Accordingly, data are analyzed in an one-day interval and the different data sets 
were added together after calibrating their a-peaks to be at the same energy point. The 
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measured pulse-height variation as a function of temperature is depicted in Figure ^ The 
PMT gain variation is typically less than 0.1%°C~^, such the origin of this pulse-height 
variation is due to change in the intrinsic light yield of the crystal. The best fit value of 
the temperature dependence is (—0.78 ± 0.17)% °C~^, which is consistent with the results 
of -0.67% °C~^ at 20 "C for GSO light yield variation previously reported |1T7| . 

Further quantitative information can be derived from the studies of the timing of the 
events, as well as from the spectral features. 

3.1 Timing Analysis 

The timing information is used to probe the fine structures of the various cascade series, 
and the results are displayed in Table ||. 

Timing-correlation analysis was applied to the ^^^^U series where the ^^^Bi-^^^Po decay 
sequence with time-correlated P-a signatures were selected: 



2i4Bi ^ 2i4pQ + ^7^ + e" + 7's (Q = 3.28 MeV ; Ti = 19.8 min) , 
2i4po _^ 2i0p^ + a (Q = 7.69 MeV ; Ti = 164 fis) . 

Since the /3-decay energy is continuous up to the end-point energy, the initial events 
were selected with energy from 0.1 MeV to 3.3 MeV, while the second was within Scte of 
the mono-energetic a-energy. An exponential function at nominal decay half-life plus a 
constant background was employed in the fit. The inclusion of a flat background was to 
account for accidental coincidences which become significant when the sought impurity 
concentration is tiny. The results are shown in Figure ^a and summarized in Table ^ 
The measured activity is (1.34 ± 0.22) mBq/kg. 

In the ^^^Th series, the ^^°Rn-^^^Po decays: 

220Rn -^ 2i6pQ + ^ (Q = g 29 MeV ; Ti = 55 s) , 
2i6po ^ 2i2p^ + a (Q = 6.79 MeV ; Ti = 0.15 s) , 

can be tagged by an a-a delayed coincidence. Pairs of adjacent events which have the 
correct electron-equivalence energies within 2aE were selected. Since the accidental rates 
in this case are comparable to the decay lifetime of the cascade sequence, the measured 
time-correlation would be faster than the nominal lifetime, and the background would be 
a slow exponential rather than a constant. The measured distribution of the time intervals 
of such pairs of events is shown in Figure ^. It was fitted to two exponential functions 
with variable lifetimes. The slow component gives a best-fit value of TsIow = 1-1 s, which is 



a measurement of the the inverse of the background rate at this energy window. The fast 
component is Tfast = 204 ms, which can be understood by Poisson statistics considerations, 
where r^^^ = r^^^^i + Bkg = t^^^^i + r^^w ■ This relation was experimentally checked to 
within an uncertainty of 10%. Folding in an over-counting factor (originated, for instance, 
from that the parent-accidental-daughter sequence would be counted as more than one 
event) of 1.2 obtained from simulations, an activity of (96.0±8.6) mBq/kg for this cascade 
sequence is derived. Accordingly, contaminations from the '^^'^Th series are expected to 
be the dominant background contribution for the crystal. 

Finally, the decay sequence of ^^^Rn-^^^Po in the ^ss-q gg^igg with an a-a tag was 
likewise processed: 

^^^Rn -^ 2i5pQ + a (Q = 6.82 MeV ; Ti = 3.96 s) , 

2i5po ^ 2iip^ + « (Q = 6.39 MeV ; Ti = 1.78 ms) . 

Events with the corrected sequence of energies were selected. The fit to the data with an 
exponential decay at nominal half-life plus a constant background is shown in Figure ^c. 
The summary of the results is given in Table |^. 

The natural abundance of ^^^U in uranium is 0.7%. However, the results based on this 
analysis shown in Table ^ implies a 235-[j^238-[j j^g^^JQ ^q ]^q 4%_ This inconsistency suggests 
that the assumption of having the decay series in secular equilibrium in the crystal is not 
valid. 



3.2 Energy Response to 7 and a 

The energy calibration to photons and electrons was achieved by the standard 7-ray 
sources (^^^Ba, ^^Na, ^^"^Cs, ^^Mn). This was fitted to a response function 



/Oe\2 ^ /^\2 /Bx2 



^C,f <o 



where the energy E is in MeV. The two terms originate from the contributions of photo- 
electron statistics and electronic noise, respectively. The calibration measurements gave 
best-fit coefficients of A=9.4±1.4% and B=1.7±0.7%, which were adopted in subsequent 
analysis. The energy resolutions are 12.4%, 8.5% and 7.2% for 0.66 MeV, 1.27 MeV and 
1.73 MeV, respectively. 

Different particles (electrons, protons, a, ...), give rise to events with different ioniza- 
tion densities dE/dx. High charge densities will lead to reduced light yield or "quenching" , 
due to recombinations and other non-radiative dissipation. Accordingly, heavily ionizing 

6 



particles like a's will give less light than electrons at the same kinetic energy deposited 
at the scintillator. 

The response of the crystal scintillator to a-particles is very important to evaluate the 
intrinsic radio-impurities background. The quenching factor of the a's depends on the 
energy and can be parametrized |18 



as: 

a ■ Ea 



1 + b ■ Eq, 

where E^ and ea are the initial kinetic energy of the a and its "electron-equivalence" 
light yield, respectively. Based on the peak positions of the a-particles from ^^^Gd 
and ^''^Am (E^ = 2.14 and 5.4 MeV, respectively), the values of a = 0.19 MeV"^ and 
b = 0.009 MeV"^ were derived and adopted in the global fitting procedures discussed in 



Section |3]^. As an illustration, a-particles with 6 MeV kinetic energy, typical for decays 



in the ^ss-q g^^^ 232rjj^ series, produce an electron-equivalence light output of 1.08 MeV. 

3.3 Spectral Shape Analysis 

A detailed Monte Carlo simulation on the contributions of intrinsic radio-impurities to the 
background spectrum was performed. The simulation software was based on the EGS4 
package |]19|. The electron and photon interaction cross sections with a GSO crystal were 
created by specifying the correct composition of 2:1:5 for the Gd:Si:0 ratio and the exact 
geometry of the crystal. Electrons and photons were traced down to the energy of 10 
keV before discarding. The number of events generated is of the same range as that 
recorded in the data taking time of 1656 hours, such that the statistical uncertainties of 
the simulated and measured spectra are comparable. 

An event generator was developed incorporating all the decay schemes based on the 
standard values [|l^ in the ^^^Th and ^ss-q ggj-jgg xhe minor (< 10~^) contributions from 
the ^^^U series were not included. The sources were randomly and uniformly located in 
the crystal volume, emitting a, jS, and 7 particles at arbitrary directions. The different 
responses to photons/electrons and to a's were incorporated based on the equations in 
Section |3.2| . While a and (3 emissions deposit all their kinetic energies in the crystal, the 
self-attenuation to the 7's by the crystal are accounted for by simulations. The combined 
effects of short-lived (< 1 /xs) cascades were generated by convoluting the various decay 
channels and their energy depositions. 

A global fit was performed, adopting as free parameters the unknown activities of all 
the long-lived precursors of the ^^^Th and ^ss-q ggj^igg^ pi^g ^j^at of the a-activities from 
^^^Gd. The ^^^Th series has three families with long-lived precursors, as shown in Table ^ 



232rp]^ and ^^^Ra are unconstrained, while activities of ^^^Th and its daughters are con- 



strained by the measurements due to timing analysis given in Section p.l| . Similarly, the 
^^^U series has five families, as shown in Table ^: ^^'^U, ^^^U, ^^°Th and ^^°Pb are uncon- 
strained, while activities of ^^^Ra and its daughters are bounded within the measurement 



values in Section 3.1 



The energy spectra due to each of these precursors and their daughters were simu- 
lated at the nominal resolution and quenching parameters based on Equations [l| and ^ 



respectively, in Section |3^. The global fit was performed by varying the relative activi- 
ties between the families to match the overall spectral shape, within the energy range of 
300 keV to 1.9 MeV. Events below this range are subjected to instabilities in electronic 
noise and trigger conditions during data taking, while those with energies above 2 MeV 
are mostly due to ambient radioactivity. Statistical uncertainties for both the measured 
and simulated spectra were taken into account. 

The best fit spectrum is displayed alongside with the measured spectrum in Figure ^ 
An illustration to the contributions from individual components to the spectral shape is 
shown in Figure ^ a's and /?+7's for the (a) ^^^Th and (b) ^^^U series, respectively. The 
measured activity for ^^^Gd a-activity (2.9 Hz in this crystal) corresponds to an isotopic 
abundance of 0.29%, larger than the nominal values of 0.2%. An alternative explanation 
for the enhancement of the a-peak at the electron-equivalence energy of 0.39 MeV can 
be the presence of ^^^Sm contamination at the 8 x 10"''' g/g level. This isotope decays 
by a-emission at an energy of 2.23 MeV and a half-life of 1.06 x 10^^ yr. The two cases, 
however, are experimentally indistinguishable. 

The activities for the families within the ^sa^jj^ ^^^ 238^ series are shown in Table ^ 
and Table |], respectively. There are two ambiguities in the exact interpretations for the 
physical processes but the difference do not have big effects on the practical concern of 
accurate background subtraction: 

I. There are three "lone" a-emissions from precursors (^^^Th, ^^^U, ^^°Th) which are the 

only member within their families. Since they are not accompanied by other decays 
and their energies are equal within 1.7 ue, the constraints from the fits are weak, 
indicating that the signatures are experimentally equivalent. Only a total activity 
of 217 mBq/kg can be derived. However, assessing from the observations that the 
crystal is dominantly contaminated by the other isotopes further down the ^^^Th 
series, it can be expected that much of the activity can be attributed to the a-decays 
of ^^^Th. The component spectra in Figure ^ assumes equal sharing among the three 
isotopes for illustration purpose. 

II. Similarly, in the energy range of interest, the ^^^U family gives two /5-decays with 



end-point of 2.2 MeV and an a-decay at 4.27 MeV. These can in principle mimic 
the lone /3's from ^^^Ac and the lone a's discussed above, both of which are at the 
same energies. Again, since the isotopes further down the ^^^Th series are known 
to be the dominant sources, the entire /3-strength was assigned to the ^^^Ac decays 
in Table § and Figure ^ 

It can be seen that the radio-isotopes within the two series are not in exact secular 
equilibrium. Nevertheless, a consistent picture can be drawn from the data where the 
^^^Th series dominate the background at a 100-1000 mBq/kg level, while the ^ss-q ggj^jgg 
are much less at the <10 mBq/kg range — except for ^^''Pb, the only long-lived daughter 
isotopes from radon contamination, which is at an activity level of 200 mBq/kg. 

Events below 100 keV originate from the /3-decay of ^^^Ra and ^^°Pb. The peak at 
390 keV is due to a-decays of ^^^Gd. The various a-decays in both the ^^^Th and ^ss-y 
series populate the broad peak around 1 MeV, defining the higher and lower energy edges, 
respectively. The shoulder between 1.4 and 2.0 MeV is due to the /5 — 7 cascades from 
^^^Ac and ^^^Bi in the ^^^Th series. The higher energy events above 2.0 MeV comes from 

(a) the /5-decay with 7-emissions of ^°^T1, which has a decay Q-value of 5.0 MeV, and 

(b) the convoluted energy oi f3 + a within the 600 ns integration-time in the ^^^Bi— ^^^Po 
cascade. The intermediate ^^^Po half-life is short (299 ns) so that 75% of decay-sequence 
are combined together as a single event. 

The relative residual (that is, [Measured— Simulated] /Measured) as a function of en- 
ergy is depicted in Figure |^a. Only statistical uncertainties are shown. It can be seen that, 
from 300 keV to 1.9 MeV where the fit was performed, the simulated spectrum matches 
well to the measured background to better than the few % level, except in the "trough" 
region between the two peaks at about 600 keV. The measured spectrum between 700 keV 
and 1.9 MeV, which is dominated by the a and (3 events for both '^^'^Th and ^^^U compo- 
nents and relevant for double beta decay analysis discussed in Section |3]^, is fitted well 
by the simulated spectra. Among them, the a-dominated bump (700 keV to 1.3 MeV) 
gives larger dispersions in their residuals. This can be understood by the fact that the 
detailed spectral shape in this region is due to the convolutions of many «-lines folded 
in with their respective quenching factors. In comparison, the shape in the 1.4 MeV to 
1.9 MeV range is due to a high energy a-decay (with energy 8.78 MeV from ^^^Po) then a 
/3-decay (Q-value 2.13 MeV from ^^^Ac), leading to more accurate and robust matching. 
The big excess above 2 MeV over contributions from intrinsic radioactivity confirms that 
ambient background, due to the 7-rays from ^os^j^ ^^-g dominant in this region. It can be 
understood from that low energy ambient photons suffer larger suppression factors due 
to the Nal(Tl) anti-Compton detector. 



3.4 Double Beta Decay Analysis 

The isotope ^^°Gd is a (3(3 emitter at a Q-value of 1.73 MeV. The (3(3 hfetime hinits 
from the present GSO measurements are not as sensitive as the various other double beta 
decay experiments H, such as the best sensitivities achieved in ^^Ge with high-purity 



germanium detectors ||2^. The analysis was nevertheless performed for completeness and 
for studying how the background subtraction results can be applied in a realistic case. 
Different approaches for setting the limits are compared. 

The measured raw background rate in the present shielding and veto configuration in 
the vicinity of the end-point is 1 kg~^keV~^hr~^. From the results in Section ^, all of 
these events can be accounted for to be due to intrinsic radioactivity, mostly /3-decays, 
at an uncertainty level of better than a few %. The contributions from these known 
sources were subtracted from the raw spectrum based on 1656 hours of data, and the 
residual spectra around the Oz//3/3 end-point is shown in Figure ||b, showing no particular 
structures. Statistical error bars were adopted. A x^ per degree of freedom (x^/dof) of 
41.9/25 was obtained for a "zero residual background" hypothesis between 1.2 MeV to 
1.9 MeV. The best-fit amplitude for a possible Gaussian signal at the expected resolution 
above zero-background at the 1.73 MeV end-point is 66 ± 1346. Assuming Gaussian 
errors, the upper limits for the signal (S^) are ^l^^^ < 699(1775) at 68(90)% CL. 

It is illustrative to compare this with the other values of S^^^ derived by alternative 
methods, which are based on the assumption that the background is a smooth spectrum at 
the end-point region, without prior knowledge of the actual background components [|T3] . 
There are 671545 events observed within the ±1 a (efficiency 0.68) region around the end- 
point, //these are all due to background, the signal must be smaller than the statistical 
fluctuations such that, folding in the efficiency factor, 

Su(90% CL) < (1.28 X VCounts)/0.68 = 1543. 
Alternatively, the spectrum can be fitted to an ad hoc polynomial and the residual at the 
±1 o" region is 39 ± 1159, giving Su(90% CL) < 2239 . Therefore, it can be seen that the 
background identification and subtraction procedures adopted in this work gives similar 
performance to subtraction schemes based on unknown background source and assumed 
spectral shape. The present scheme provides further justifications and robustness to the 
background suppression, and is apphcable for signals other than the easily distinguishable 
energy peaks (such as continuous spectra) as well as at regions where the background 
spectra do exhibit certain structures (indeed, there is a bump at 1.6 MeV due to a high 
energy a-decay). The fact that the present uncertainty is comparable to the statistical 
fiuctuations (vCounts) indicates that, up to a suppression factor of 10"^, the procedure 
is still statistics-limited and there are room to enhance the sensitivities further with more 
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data and longer simulation time. 

The half-life limit can be expressed as 

T- > .og.2 X '- -IS""'" . (3) 

2 bu 

Putting the experimental values of the target number N(^^°Gd) = 1.07 x 10^^, and the 
data taking time of 1656 hours(0.19 yr) and the values of Su derived above, we obtain 

T°"'^^ > 2.0(0.8) X 10^0 y at 68(90)% C.L. , 

2 

consistent with and at a comparable range to the previous results: Ti > 3.0 x 10^° y 
at 68% C.L. ig and Ti > 8.2 x 10^° y at 90% C.L. [|14|. These limits apply to both 



the 0^ -^ 0^ ground state to ground state transition as well as the O"*" — > 2"*" (87 keV) 
transition to the excited state, since the 87 keV 7-rays emitted are fully absorbed by 
the crystal. The corresponding limits of the effective neutrino mass parameter, based on 
calculations in Ref. [H, are (m^) < 65(103) eV at 68(90)% CL. 

Limits can also be set for the two-neutrino double beta decay (2i//3/3) and the neutrino- 
less double beta decay with Majoron emission [xPP) channels. We took the conservative 
approach of deriving the limits only with the residual data in Figure Pd between the 
1.2 MeV and 1.9 MeV region. The x^/dof is good for a zero-residual hypothesis using sta- 
tistical error bars alone, indicating that the background in this region is well-understood. 
The efficiency factors for the 2z//9/9 and x/3P channels to be within this energy range are 
2% and 46%, respectively. 

The best-fit of a 2i'f3f3 spectra over zero-residual in this energy range gives the optimal 

activity of (34620 ± 67167) events, from which S^^^^ < 66188(120593) at 68(90)% CL, 

leading to 

Tf ^^ > 2.2(1.2) X 10^^ y at 68(90)% C.L. 
2 

Similarly, a fit to the xPP spectra gives the best value of (— 1092±4138) events. Restricting 
a possible signal to be positive definite, the upper limits S^'^^ < 4676(7241) at 68(90)% CL, 
is derived, giving rise to 

Tf^ > 3.0(1.9) X 10^^ y at 68(90)% C.L. 
2 

Though these 2i//?/? and xPP limits are only modest compared to those of the other double 
beta decay candidates, they represent nevertheless more than an order of magnitude 
improvement over previous published results on ^^°Gd |T^. The improved sensitivities 
can be attributed to that the background must first be quantitatively understood and the 
detailed spectral structures accounted for to allow them be subtracted off from signals 
which are themselves continuous spectra without specific structures. 
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4 Discussions and Summary 

The analysis methods and results discussed in this article give indications of the potential- 
ities and problems of GSO crystals in neutrino experiments. They also provide training 
grounds in developing the schemes for the background understanding and suppression for 
similar class of experiments. 

Low background experimentation techniques must be used to maximize the sensitiv- 
ities in a realistic experiment. To get to the relevant 10^^ y sensitivity range for double 
beta decay searches, one needs 10^ improvement in the [(target mass x time) /background] 
value. With a realistically achieve-able target mass of 0(10 ton) and running time of 
0(1 year), the background has to be further reduced by at least factor of 10^, via more 
elaborate shielding configurations and data taking at an underground laboratory, as well 
as through careful choice of raw materials with reduced ^^^Th and ^ss-q activities, es- 
pecially the contaminations due to ^^^Ra. Measurements with raw background rate of 
0.04 kg~^keV~^hr~^ at the Oz//5/5 end-point were reported |0, such that factor of hun- 



dreds in improvements have still to be worked on. An alternative or perhaps better can- 
didate for double beta decay experiments with crystal scintillators ^^^Cd with cadmium 
tungstate (CdW04) crystals |2|. 

Solar neutrino experiments with ^^°Gd typically involve an event rate of 0(0.1 ton~^day~^) 
and hence imply more stringent intrinsic radio-purity level requirements. The specific 
delayed-correlation signatures discussed in Section |^ provides powerful means to suppress 
background due to accidental coincidence. Design studies ||^ indicate that the purity 
levels of < lO^^*' g/g are required for ^ss-q g^j^^ 232rj]-^ series and 10~^^ g/g for ^ss-q (^^j^ere 
/3-decays of ^'^^Th with a delayed 7-emissions may mimic the neutrino signatures). The 
results in Table ^ show that improvements to the 10^ and 10 levels relative to the activi- 
ties measured in this work for the ^^^Th and ^^^U series, respectively, are necessary. This 
timing analysis results in Table |^, however, only gives the contamination levels of ^^^Ac 
and its short-lived daughters within the ^ss-q gg^-igg i^^^ does not provide direct measure- 
ment of the ^^^U concentration. If secular equilibrium is assumed, an improvement of 
better than 10^ is necessary. The electron-equivalence light output of 392 keV for ^^^Gd 
a-decay is in between the expected range of t'e(pp) and z/e(^Be). Effects of the tails of the 
energy distribution to the < 10~^ level have to be studied. There are isotope-separation 
schemes to deplete ^^^Gd [^, thereby reducing this potential background. In addition. 



care must be exercised in selecting the sources of cerium the dopant materials in order to 
avoid contaminations from its long-life isotopes: ^^^Ce (Ti = 140 days) which decays by 
electron capture emitting 7-rays of 166 keV, and ^^^Ce (Ti = 284 days) which decays by 
/3-emissions at a Q- value of 320 keV. 
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A GSO scintillating crystal detector which can satisfy the requirements for solar neu- 
trino experiments will also be a powerful device for double beta decay searches that can 
operate simultaneously. 

One of the problems for the construction of big detector systems based on GSO crystals 
has been the cost. There are ongoing efforts from the industries to produce and market 
this crystal at a similar price range as BGO [^. If achieved, this crystal will be an 
attraction option for electromagnetic calorimeter in a high-rate environment as well as 
for medical imaging applications. 

An independent motivation of this work is to study the methods for the measure- 
ments and subtraction of background due to intrinsic radio-isotopes. In a scintillating 
crystal detector with minimal passive volume and where a three-dimensional fiducial vol- 
ume definition can be achieved, the ambient background in the sub-MeV range will be 
highly reduced and the primary experimental focus will be one on the understanding and 
suppression of internal background [Q]. In addition, forthcoming big liquid scintillator 
experiments like BOREXINO [^ and KamLAND [^ for low energy solar neutrino and 



long-baseline reactor neutrino experiments will need to identify and measure the back- 
ground processes quantitatively to justify their signals, in a similar spirit to this work. 
The GSO crystal used in this measurement, having relatively high contamination levels, 
provides a convenient platform to develop the analysis procedures and to understand the 
various systematics involved. Crystals with much better radio-purity level are available. 
For instance, the secular equilibrium levels for the ^^^Th and ^^^U series in CsI(Tl) are 
below the 10"^^ g/g range |^. 

It has been shown that a subtraction factor of 10~^ is possible, at least in the case 
where the region is dominated by a smooth /?-decay background. The background of 
232rp]^ and ^^^U in a GSO crystal treated in this work are in fact more difficult to han- 
dle from the background subtraction point of view. There are multiple a-decays to be 
convoluted together with their quenching factors to be taken into account, detection ef- 
ficiency factors for the 7-emissions which depends on the details in detector geometry, 
as well as convolutions of various signals due to short life-time cascades. Other simpler 
decay schemes are expected to be easier to handle. For instance, the background due 
to /3-decays of ^°K can be more accurately accounted for, since these events are uncor- 
related with the other processes such that the background spectral shape will be simply 
the well-known /3-spectra. Moreover, the activity can be independently obtained, besides 
being derived from the best-fit value, through the measurements of the related 1.46 MeV 
7-decays of ^°K. Furthermore, crystals like Nal(Tl) and CsI(Tl) Q, which have excellent 
pulse shape discrimination properties for 7/a separation and where the a-particles suffer 
less severe quenching, can provide much better measurements and identifications for the 
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a-decays from the various channels, providing further powerful diagnostic tools to probe 
the structures within the various cascade series. 

The authors are grateful to T. Murakami for the loan of the GSO crystal, and to C. 
Rangacharyulu for discussions on the various decay schemes of radio-isotopes. The work 
is supported by contracts NSC88-2112-M-001-007 and NSC89-2112-M-001-028 from the 
National Science Council, Taiwan. 
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Properties 


Nal(Tl) 


CsI(Tl) 


CaF2(Eu) 


BGO 


GSO 


Density 


3.67 


4.51 


3.19 


7.13 


6.7 


Relative Light Yield 


i.oot 


0.45 


0.50 


0.15 


0.2 


Emission Peak Am (nm) 


415 


530 


435 


480 


440 


Decay Time (ns) 


230 


1000 


940 


300 


30-60 


Refractive index at Am 


1.85 


1.80 


1.44 


2.15 


1.85 


Hygroscopic 


yes 


slightly 


no 


no 


no 


Melting point (°C) 


651 


621 


1418 


1050 


1900 



t Typical light yield for Nal(Tl) is about 40000 photons per MeV. 

Table 1: Characteristic properties of several common crystal scintillators and their com- 
parison with the GSO crystals 



17 



Decay series precursor 


232Th 


238U 


235U 


Actual decay sequence 


220Rn-2i6po 


214Bi-214po 


2i9Rn-2i5po 


Signatures 


a — a 


(5-a 


a — a 


Bounds on time interval (ms) 


[10, 1000] 


[0.13, 1.3] 


[0.32, 4.8] 


Half-life (ms) 


150 


0.164 


1.78 


Duration of Data Taking (hours) 


48 


559 


559 


No. of events 


28259 ± 2543 


2651 ±445 


1015 ±588 


Efficiency of Cuts 


0.82 


0.57 


0.78 


Over-counting factors 


1.2 


1.0 


1.0 


Activity (mBq/kg) 


96.0 ±8.6 


1.34 ±0.22 


0.37 ±0.21 


Concentration of Precursor (g/g)^ 


(2.9 ±0.3) 


(1.1 ±0.2) 


(4.4 ±2.5) 




xlO-^ 


xlO-io 


XlO-12 



^ Assuming secular equilibrium 

Table 2: Summary of results in the determination of the concentration of the thorium and 
uranium series in the GSO crystal based on energy-timing correlations. The statistical 
uncertainties are quoted. 
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Family in 


Isotope 


Half-Life 


Decay Mode 


Q- value 


Activity 


232 Th series 








(MeV) 


(mBq/kg) 


232Th 


232Th 


1.405 X 10^° y 


a 


4.013 


t 


228Ra 


228Ra 


5.75 y 


P 


0.046 


733 t 




228AC 


6.15 h 


(3 


2.127 




228Th 


228Th 


1.913 y 


a 


5.423 


107 




224Ra 


3.66 d 


a 


5.685 






220Rn 


55.6 s 


a 


6.288 






216Po 


0.145 s 


a 


6.778 






212pb 


10.64 h 


13 


0.574 






212Bi 


60.55 min 


a (64%) 
13 (36%) 


6.090 

2.254 






212po 


0.299 fis 


a 


8.784 






208 fjj 


3.053 min 


P 


5.001 




Family in 


Isotope 


Half-Life 


Decay Mode 


Q- value 


Activity 


23^ U series 








(MeV) 


(mBq/kg) 


238U 


238U 


4.47 X 10^ y 


a 


4.27 


<2t 




234Th 


24.1 d 


13 


0.183 






234pa* 


1.18 min 


P (99.87%) 
IT (0.13%) 


2.207 
0.08 






234pa 


6.70 h 


P 


2.287 




234U 


234U 


2.45 X 10^ y 


a 


4.856 


t 


230Th 


230Th 


8.0 X 10^ y 


a 


4.771 


t 


226Ra 


226Ra 


1600 y 


a 


4.871 


< 2 




222Rn 


3.83 d 


a 


5.591 






218Po 


3.05 min 


a 


6.115 






214pb 


26.8 min 


(3 


1.024 






214Bi 


19.7 min 


f3 


3.270 






214po 


163.7 /is 


a 


7.834 




210pb 


210pb 


22.26 y 


P 


0.063 


198 




210Bi 


5.01 d 


P 


1.161 






210Po 


138.4 d 


a 


5.408 





^ Experimentally Equivalent. Total 217 mBq/kg. See Text. 

■'■ Experimentally Equivalent. Assign weight to 2^8 j^a. See Text. 



Table 3: Activities of the radio-isotopes in the 232^]-^ and 238^ series, based on best-fit 
parameters to the background spectrum. 
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Figure 1: Schematic drawing of (a) the experimental setup showing an active Nal(Tl) 
veto scintillator around the GSO crystal, together with a plastic scintillator for cosmic-ray 
veto on top and 4 tt lead shielding with at least 10 cm in thickness ; (b) the electronics 
and data acquisition system. 
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Figure 2: Comparison between the measured spectra from 1656 hours of data with the 
simulated spectra based on the best-fit rates from the various components. 
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Figure 3: The measured light yield variation of the GSO crystal, using the ^^^Gd a-peak 
as reference, as a function of temperature during the course of 1656 hours of data taking. 
The best-fit dependence is overlaid as the straight line. 
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Figure 4: Distribution of the time intervals between (a) the f3 decay of ^^^Bi and the 
following a decay of ^^^Po, from 559 hours of data; (b) the a decay of ^^°Rn and the 
subsequent a decay of ^^^Po, from 48 hours of data; (c) the a decay of ^^^Rn and the 
following a decay of ^^^Po, from 559 hours of data. The solid lines shows the fits of the 
data as discussed in the text. 
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Figure 5: Simulated energy spectra due to the a and /? + 7 contributions from the (a) 
^^^Th and (b) ^^^U series. 
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Figure 6: (a) Ratio of the residual (measured— simulated) to the measured background 
spectrum for 1656 hours of data, (b) Residual events at the energy range relevant for 
double beta decay analysis. Statistical error bars are shown. 
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